Free surface flows with large slopes: beyond lubrication theory 
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The description of free surface flows can often be simplified to thin film (or lubrication) equations, 
when the slopes of the liquid-gas interface are small. Here we present a long wavelength theory 
that remains fully quantitative for steep interface slopes, by expanding about Stokes flow in a 
wedge. For small capillary numbers, the variations of the interface slope are slow and can be 
treated perturbatively. This geometry occurs naturally for flows with contact lines: we quantify the 
difference with ordinary lubrication theory through a numerical example and analytically recover 
the full Cox-Voinov asymptotic solution. 
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Free surface flows are encountered in many everyday 
life and industrial situations, ranging from soap films, 
sliding drops to paints and coatings [lj, |2| • The hydrody- 
namic description of these 'free boundary problems' still 
provides a challenge of great fundamental and techno- 
logical interest. The difficulty lies in the intricate cou- 
pling between the liquid-gas interface and the flow inside 
the film, which gives rise to a broad variety of insta- 
bilities and interface morphologies 

HiaSlH In 

the case of 'thin' films, for which horizontal and vertical 
length scales are well separated (Fig. QJi), the problem 
is greatly reduced through a long wavelength expansion 
0,H, . At low Reynolds numbers this so-called lubrica- 
tion approximation yields a single nonlinear equation for 
the evolution of the interface profile h(x, y, t), and forms 
the accepted theoretical framework for free surface flows. 
This reduction is possible whenever surface tension (7) 
dominates over viscosity (77), so that the capillary num- 
ber Ca = r\U* /7 serves as a small parameter; U* denotes 
the velocity scale of the problem. 

The standard formulation of lubrication theory, how- 
ever, has a severe drawback: it is only valid for small 
interface slopes, i.e. |V/i| <C 1. While it is generally 
believed that the lubrication equation yields good qual- 
itative predictions for larger slopes as well, one has to 
be careful with quantitative comparisons. This is par- 
ticularly important for the problem of moving contact 
lines for which viscous forces tend to diverge as h — ► 
[ToL HH . The microscopic mechanisms that release this 
singularity are highly disputed 0, 0, Q, 0, 0] , and 
call for a fully quantitive description of experiments that 
often involve large contact angles. As it is practically 
infeasible to resolve the full hydrodynamic problem on 
all relevant length scales, ranging from molecular to mil- 
limetric, a simplified theory for finite slopes would be 
extremely valuable. 

In this paper we present a generalization of the lubrica- 
tion theory for free surface flows at low Reynolds numbers 
that remains exact for large slopes. The crucial observa- 
tion is that in the limit of small Ca, capillary driving 
requires slow variations of the interface curvature, but 
there are no restrictions to the steepness of the interface. 



(a) 





z=h(x,v) 


gas 




y 


liquid 



h* 



(b) 




FIG. 1: (a) The usual lubrication approximation is valid 
whenever the liquid-gas interface is slowly varying along the 
horizontal coordinates, and is thus restricted to small slopes, 
(b) Considering wedge-like profiles with slowly varying slope 
9(x), we derive the long wavelength theory for steep slopes. 



We therefore consider the flow in a wedge with a finite, 
but slowly varying opening angle 8 - see Fig. This 
geometry naturally occurs for contact lines. Treating the 
variation as a perturbation around a straight wedge then 
yields the equation for the interface profile h(x, y, i). 

This theory, summarized by Eqs. i|14ll5|) . remains fully 
quantitative for large slopes when the curvature k<1 /h, 
while it enjoys the same mathematical structure as the 
usual lubrication theory. We furthermore show that the 
equation reproduces the asymptotics for 9{x) as com- 
puted by Voinov 01 an d Cox , in a relatively straight- 
forward manner. However, the present work reaches be- 
yond asymptotic relations: it describes all intermediate 
length scales as well, and allows to incorporate other 
forces such as a disjoining pressure or gravity. 

Lubrication theory for free surface flows - Before ad- 
dressing the problem of finite slopes, let us first briefly 
revisit the lubrication approximation. In the limit of zero 
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Reynolds number, the flow of incompressible Newtonian 
liquids is described by Stokes equations 

V-v = 0, (1) 
-Vp + ?7Av-V$ = V x Av = 0, (2) 

where p(x, y, z) and v(a;, y, z) represent the pressure and 
velocity field respectively while we consider body forces 
that derive from a potential <I>. (Time dependent profiles 
will be discussed below.) The equation is complemented 
with the boundary condition of Laplace pressure at the 
free surface 

p(z = h) =p --yK, (3) 

where k is the interface curvature. This gives rise to 
an intricate nonlinear coupling between the shape of the 
interface and the flow inside the film, which has to be 
resolved self-consistently. We limit the discussion to the 
case where the gas is hydrodynamically passive and we 
take a zero shear stress condition at the free interface. 

In the limit where Ca <C 1, the free boundary problem 
can be greatly reduced through the lubrication approxi- 
mation: surface tension is sufficiently strong to drive the 
viscous flow through only minor variations of the shape 
of the free surface (Fig. GJi). One thus expects that the 
interface profile and the velocity field are slow functions 
of the horizontal coordinates, so that 




while we write 

i \ ~ ( x y z \ 

= vo + e vi + e 2 v 2 + ■ ■ • . (5) 

All lengths have been rescaled by the typical film thick- 
ness h* , and e = Ca 1 ^ 3 is the small parameter of the 
expansion. The strategy is to solve Eqs. lHIMI) order by 
order in e. Here we briefly sketch the approach; for a 
more detailed derivation we refer to 0- 

If we let e go to zero, the departure of h(x,y) from a 
horizontal interface becomes increasingly small. Hence, 
the velocity profile converges towards the parabolic 
(Poiseuille-like) profile in this limit, so that 

* -£(,-(,_£)■). m 

Formally, this dominant flow can be obtained from 
Eq. J2J a t order e°. Note that Vo still evolves on a long 
scale through its dependence on h(x,y). 



In deriving Eq. © we used the boundary conditions 
of no-slip at z — 0, and zero shear stress at z = h. The 
prefactor has been chosen such that U = l/hf Q dzvo 
represents the depth-averaged velocity in the frame at- 
tached to the plate. 

Since the dominant viscous forces in Eq. J2J) arise from 
Vo, we do not need to solve for the higher order veloc- 
ities to obtain the equation for h(x,y). At leading or- 
der, Eqs. J2J an d © reduce to the celebrated lubrication 
equation 

U/U* 1 

VAft = 3Ca + -V$U=fc- ( 7 ) 
Combined with the depth-averaged continuity equation, 

d t h + V ■ (hV) = 0, (8) 

it provides the common theoretical framework for free 
surface flows, both in the scientific community as well as 
for industrial purposes. The lowest order terms that are 
neglected are of order e 2 , so the expansion is valid in the 
limit of small slopes, |V/i| 2 <C 1. 

Theory for large slopes - We now perform a simi- 
lar long wavelength expansion for wedge-like geometries, 
such as depicted in Fig.^x The crucial physical ingredi- 
ent underlying the expansion is that surface tension dom- 
inates over viscosity, i.e. Ca < 1, so that variations of 
interface curvature are slow with respect to the relevent 
length scales. In principle there is no restriction to the 
slope of the interface: the only requirement is that the 
opening angle 8 is slowly varying (Fig-UJa). We therefore 
consider profiles with 

e(x) = e(e£), (9) 

and perform an expansion in e. Note that 8 itself is of or- 
der unity. Here we introduced the length scale x*, which 
is a typical distance to the 'origin' of the wedge - we show 
below that the analysis remains self-consistent as long as 
hd x 8 -C 1. For simplicity we discuss two-dimensional 
profiles so we omit the y-dependence. 

In the spirit of the lubrication approximation, we ex- 
pand the velocity as 

v(x, z) = v + e vi + e 2 v 2 + • • • , (10) 

and solve for the dominant flow vo . Again, vo is obtained 
in the limit e — ► 0, which in this case corresponds to a 
wedge of constant opening angle. Hence, the problem for 
Vo reduces to Stokes flow inside a straight wedge, which 
is easily solved analytically 0, 0] . 

To make this more explicit, we introduce local cylin- 
drical coordinates, r(x, z) and 4>{x, z), which are defined 
by the locally tangent wedge (Fig. |2K) 
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FIG. 2: (a) Definition of the cylindrical coordinates r(x, z) 
and <p(x, z) in the locally tangent wedge of angle Q(x). (b) The 
basic velocity vo corresponds to flow in a wedge of constant 9, 
sketched in the frame comoving with the interface, Eq. 112t . 
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The ^-dependence appears through h(x) and 9{x). Writ- 
ing the velocity as a function of these coordinates 



v(x,z) = v(r(x,z) , (j>{x,z)) 



(11) 



and expanding Eqs. (11121) . one indeed finds that the order 
e° reduces to the problem of a straight wedge: variations 
of 9 show up at higher orders. For the dominant flow we 
can thus use the results of [ToL \v\ , 



(V )r = U 
(V )* = U 



(cos <f> — (j) sin <j>) sin 9 — 9 cos 9 cos < 

6> — cos 9 sin 6* 
6* sin <b cos — 6 cos (A sin 



(12) 



which hold in the frame comoving with the interface. 
This flow has been sketched in Fig. [2t>- Here we used 
the conditions of a vanishing shear stress at 4> = and 
a no-slip condition at the plate, v r = —U. The latter 
condition ensures that U represents the depth-averaged 
velocity in the frame attached to the plate. 

Evaluating Av at the free surface (</> = 0), Eq. (0 
provides the leading order pressure gradients along the 
interface 



Combined with the Laplace pressure condition © this 
yields the generalized lubrication equation: 



U/U* 1 
3Ca^—F(9) + -d x $\ z=h , 
n z 7 



where 



F{9) = 



2 tan 9 shr 9 

3 9 — cos 9 sin t 



(14) 



(15) 



Since d x n ~ (e/x*) 2 , the expansion is self-consistent 

1 /9 

when taking e = Ca ' . 

Comparing this result to the lubrication equation 0, 
one observes two differences. First, the function F(9) can 
be seen as a correction factor for the viscous term: in- 
deed, one recovers F{9) — » 1 for small slopes. Second, the 
left hand side of Eq. 114JI now involves the full curvature 



{l + d x h 2 ) 3 / 2 ' 



(16) 



The expansion thus provides an equation for the inter- 
face profile that has the same mathematical structure as 
the usual lubrication approximation, but which remains 
exact for large slopes. The lowest order terms that are 
neglected are of order e = 0(x*d x 9), so the expansion 
is valid when hd x 9 <C 1. This does not mean that the 
description is limited to small heights: the height only 
appears as a reference scale to quantify the dimension- 
less curvature. 

Asymptotics for 9(x) - To illustrate the strength of 
the approach, we now show that Eq. (|14|) correctly re- 
produces the nontrivial asymptotic solution for 9(x), as 
obtained by Voinov and Cox 0] for advancing con- 
tact lines. Anticipating the well known result, we express 
the relation between 9 and x as 



g(9) = Caln(a:/a;o) 



(17) 



and solve for g(9) using Eq. I|14|l with $ = 0. Differenti- 
ating g{9) with respect to x one finds d x 9 = Ca/(xg'), so 
that the curvature can be written as k = Ca cos9/(x g'). 
Combining this with the expansion h(x) = xia,\\9(x) + 
O(Ca), Eq. (|14J) becomes to lowest order 



Ca cos 9 
x 2 g' 



3Ca 



c 2 tan 2 9 



F(9). 



(18) 



Here we used the fact that advancing contact lines move 
along the negative x-direction in the frame attached to 
the plate, so that J7 ac i v . = — U* . Eliminating g' , we read- 
ily recover the famous result of 



drP\<t>=0 = - 



suit 



(13) 



9(0) 



du ■ 



■ cos u sinu 
2 sin u 



(19) 
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FIG. 3: Numerical solution of a flat film dragged downwards 
by gravity in the presence of a contact line. The film thickness 
ftflim is uniquely determined by imposing a contact angle 6o 
at a microscopic scale, here taken ho = 10~ 5 (lengths are 
expressed in the capillary length y^fjpg)- At large angles, 
the results of Eq. 1141 (solid) provide significant corrections 
with respect to lubrication theory Eq. (7J (dashed). 



Discussion - We have derived the long wavelength 
expansion for free surface flows in the case of steep in- 
terface profiles. This provides a significant improve- 
ment with respect to the usual lubrication approxima- 
tion, whose validity is restricted to small slopes. The 
resulting theory has the same mathematical structure as 
the lubrication equation and is thus easily adapted to 
existing codes and methods. The most natural applica- 
tion of our work is found in wetting flows that involve 
large contact angles. This is illustrated in Fig. |31 where 
we computed the thickness of a flat film flowing down a 
vertical plate in the presence of a receding contact line. 
The effect of gravity is accounted for through $ = —pgx, 
which introduces the capillary length scale l 7 — y/j/pg- 
For a given thickness hfu m = Z 7 V3Ca, there is a unique 
solution that connects to the contact line [l^. By nu- 
merically solving the interface profile down to a micro- 
scopic (molecular) height ho, we can thus identify the 
slope v ery close to the contact line, denoted by the an- 
gle #o liil- This yields a unique relation between the 
slope imposed at a microscopic distance from the con- 
tact line and the macroscopic film thickness. The results 
obtained from numerical integration of Eq. I|14|) (solid) 
displays significant quantitative differences with the pre- 
diction by lubrication theory (dashed) at angles > 30°. 

Although in principle the expansion is not limited to 
contact lines, steep slopes are typically attained through 
curvatures hh" ~ 1 in flows without contact lines. This 
lies beyond the strict validity of the expansion and one 
expects corrections due to curvature of the interface. 
These corrections can in principle be treated perturba- 
tively as well, since Ca <C 1 implies weak variations of 
curvature. 
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